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waves and generate infra-red divergences. In this paper we will show that 
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fields which do create in and out Fock states. In a companion paper we 
demonstrate that the Green's functions of these fields have a good pole 
structure describing particle propagation. 
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1 Introduction 



The successes of hadronic physics are to a great extent based on our use of partonic 
variables in the ultraviolet (UV) domain. However, in the realm of soft physics, where the 
dynamics of quarks and gluons underlie jet formation and hadronisation, we simply do 
not know what variables to use. Strong interaction phenomenology could be a primary 
source of information here and this is, to a large extent, based upon effective quark 
models. However, the 'quarks' which compose these models are, although undoubtedly 
useful in very many ways, in no sense organic to Quantum Chromo dynamics (QCD). 
This remains the case when such models are further mutated in the hope of making them 
more closely resemble the low energy dynamics of QCD, e.g., by introducing 'gluons' 
into effective quark models. Whatever these effective variables are, they are not simply 
the Lagrangian fields of QCD. 

Knowledge of the correct infra-red (IR) variables in QCD would help us to obtain 
a deeper understanding of the structure of the underlying field theories of the standard 
model, to obtain some feeling for the validity and the (a priori unknown) limits of 
phenomenological models in particle physics and, most importantly, help us understand 
hadronisation and confinement. This paper is devoted to the construction of physical 
charged fields using only the gauge theories of the fundamental interactions. Faced 
with such a major task, we cannot hope to solve it immediately and our procedure will 
be essentially perturbative. Our aim is to develop methods and insight which may be 
carried over and extended into the non-perturbative domain. We intend to show below 
that this aim is not completely unrealistic. 

Let us first recall how the problems associated with the use of unphysical variables 
show themselves in the standard description of scattering. In the LSZ approach it is 
assumed that at times long before or after any scattering process the fields entering or 
emerging from the vertex do not interact with each other any more. The experimental 
fact of hadronisation tells us that this cannot be true if our asymptotic fields are just 
taken to be the Lagrangian quark fields. Although we can describe the hadronisation 
process in one of the many models, this is necessarily matched on to partons emerging 
from the scattering process. This has even been elevated to the status of the principle 
of local parton hadron duality |lj>0- However, this matching onto free partons has a 
price: quarks and gluons still interact with each other even when widely separated. 
The various IR divergences in on-shell Greens functions and in S-matrix elements are 
a direct consequence of the neglect of such interactions. In fact even in Quantum 
Electrodynamics (QED), where there is a 1/r fall-off in the interaction between the 
charges which we will generically refer to as electrons, it has long been known p|-|T0 



that this fall-off is still too slow to permit us to ignore the interaction and its neglect is 
responsible for the IR catastrophe. 

There are various responses to this situation: the one most widely, if tacitly, employed 
has been to give up any description of on-shell Greens functions and S-matrix elements 
and to restrict oneself to IR safe quantities, such as, in QED, cross-sections summed 
over all final state soft photons. This response is at once theoretically radical (discarding 
the S'-matrix) and practically conservative (we should only attempt to talk about this 
class of measurable quantities). It may be noted that this posture does not address 
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the questions initially raised above, viz. what are the correct IR degrees of freedom? A 
second approach is to redefine the division of the Hamiltonian into free and interaction 
terms so that the new interaction can indeed be neglected at asymptotic times. This 
approach has been studied for QED by various authors [0,11,0,0 but has not been 
successfully extended to QCD. As stressed by the authors of ||, taking this route entails 
the loss of a particle description already at the level of QED! (The perturbative signal 
for this is that the on-shell QED Green's functions so obtained, although IR finite, do 
not have a pole structure describing particle propagation, rather the external lines for 
electrons and other charged fields have branch functions.) A third reaction, which we 
pursue below, is to argue that the experimentally observed particles entering and leaving 
scattering processes do not directly correspond to the variables which we have in the 
Lagrangian. E.g., when an electron emerges from a scattering process it is accompanied 
by an electromagnetic field which is not accounted for by the free Lagrangian fermion, 
■0. We will refer to the inclusion of such fields as dressing the Lagrangian matter field. 
Such ideas have been investigated by many authors, see, e.g., |13|, [L4|, p!^-p8 . 



We will argue below that this last approach, choosing the correct asymptotic fields, 
indeed solves the problem in a way which is compatible with having the expected pole 
structure, retaining the S-matrix and having a particle description. Here of course we 
need to understand what are the 'correct' fields. We will present below two requirements 
for the construction of such fields and further show how the descriptions which satisfy 
these conditions may be obtained in a systematic fashion which can be extended to other 
theories and which will permit the eventual inclusion of non-perturbative dynamics. 

We have seen that the interaction in gauge theories cannot be neglected at large 
times. Now the Lagrangian fields in QED transform under gauge transformations as 

A^x) ^ A^x) + d^x) } and ^(i)^e l *VW, (1) 

so, since the coupling does not vanish, the matter field, ip, is not gauge invariant at 
large times and cannot be identified with a physical particle such as an electron even in 
the asymptotic region. Can we then construct a gauge invariant charged field at all? 

To make this more precise, we can use the tools of constrained dynamics (see, e.g., 



[p9| ). The prime requirement, as far as we are concerned in this paper, is that a physical 



field must satisfy Gauss' law, which reads 

&Fn = -e J . (2) 

Since Gauss' law generates gauge transformations, this translates into the requirement 
of gauge invariance. There is a great deal of work in the literature on gauge invariant 
variables or simply variables with simplified gauge transformation properties, see, e.g., 



30|-|32|j. However, the identification of such variables with charged particles is especially 
subtle. The main point we wish to note here is that Gauss' law tells us that charged 
matter fields cannot be separated from the electromagnetic cloud which surrounds them. 
This implies then that any description of a physical charged particle cannot be local 
since the total charge can always be written as a surface integral at spatial infinity . 
We should also point out that, since the charge density generates global (rigid) gauge 
transformations, charged particles are not invariant under such gauge transformations. 



2 



The first construction of charged fields that we are aware of was performed by 
Dirac fl3||, who noted the existence of a set of composite fields 



i{> f (x) = exp (-ie J d A zf»{x - z)A^zfj ^(x) , (3) 

which, Dirac argued, are locally (but not globally) gauge invariant for any / M so long as 
duf^w) = 6<V(w) holds. 

It is easy to see that there are many gauge invariant descriptions (many functions 
/ which fulfil this property). One of the main ideas behind the work reported here is 
that gauge invariance is not enough: a physical field must single itself out by further 
requirements. It is easy to see that some constructions are gauge invariant but highly 
unphysical: the most obvious example being the stringy ansatz 



X 



^t(x) = exp -ie / dziAi(x , z) ip(x) , (4) 



r 



which, it can be shown corresponds to an infinitely excited state where the electric 
flux is confined along the path of the string, T. The reader can find an animation 
showing the decay of such a state at http://www.ifae.es/~roy/qed.html. 

Dirac suggested the requirement that, out of the general set of functions given by 
Eq. |3|, the physical description should have the correct electric field. Concretely, he 
noted that a physical charged particle is accompanied by an electric field and argued 
that an electron should be described by 

ip D (x) = exp (-ie^fj ip(x) , (5) 

which is easily seen to be a special case of (|j). His argument for this choice was that 
the state, i/jd(x)\0), has a Coulombic electric field 

^, aj )^(y)|0>=--^ 1 ^l^ i ,(»)|0>. (6) 

47r \x — y\ 6 

To obtain this result^, he used the equal-time commutator: [Ei(x), Aj(y)] = i5ij5(x — y). 
We note that ipo is locally gauge invariant and that the factor of 1/V 2 in (|5]) shows the 
non-locality of the electromagnetic cloud around the charge, since 



4-S(x) = ~\^ XT i^L.. (7) 



\ z 4tt J \x — y\ 

Since no spatial direction is singled out, we might naively expect this to probably de- 
scribe a static charge and this is in agreement with (||). 

A natural question now is how we should interpret this gauge invariant field? Dirac 
argued that it corresponds to a charge at y together with its associated Coulombic 
electric field. Since the Coulomb field is that of a static charge, his interpretation implies 
that we know both the position and the velocity of the charge! This, as is well known 

4 More generally, e.g., in Coulomb gauge where the exponential factor reduces to unity, one needs to 
employ Dirac brackets. 
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from heavy quark effective theory (HQET), is only compatible with the uncertainty 
principle in the infinite mass limit. This will be clarified in Sect. 3 for charges with 
finite masses. 

An obvious test of Dirac's construction is to study whether these fields remove the IR 
singularities associated with on-shell static charges. The simplest test is the two-point 
function of charges described by @ in an on-shell renormalisation scheme. It was shown 
in [53| that this is indeed IR finite so long as the static point on the mass shell is chosen, 
i.e., p = (m, 0, 0, 0). (The generalisation of this to the propagator at an arbitrary point 
on the mass shell was shown in MM for scalar QED and in |35j for the fermionic theory.) 



We note that it was further shown in |36| how to extend this description of dressed 
electrons to construct gauge invariant dressed quarks in non-abelian gauge theories at 
any order in perturbation theory (see the appendix of |24|). Essentially there is a 



minimal perturbative extension of the QED result which retains gauge invariance. We 
note that Haller and his collaborators [0,0 have also obtained perturbative formulae 
which agree with this minimal extension of the static charge (|5|). 

Three important implications of this work for QCD are that colour charges are 
only well defined []39| for locally gauge invariant fields such as our dressed quarks and 



gluons; the identification of a topological obstruction to the construction of an isolated 
quark [24] and the perturbative identification [|G] of the dominant gluonic configuration 



responsible for asymptotic freedom. 

In the rest of this paper we will present a method to construct charged particles in 
QED. We will systematically obtain explicit solutions and further demonstrate that the 
charged fields we construct, since they are surrounded by the correct electromagnetic 
fields, obey a free asymptotic dynamics. It will become clear that our results reflect 
previous work on the structure of the asymptotic dynamics which is found when it is 
assumed that the fields entering and leaving scattering processes are just the free matter 
fields, however, in our approach we will retain a particle interpretation. 



In a companion paper ||41|] , II, the results which we obtain below, and the interpre- 
tation which we give to them, will be further verified. In particular, we will use our 
physical variables and study the IR behaviour and UV properties of various Green's 
functions and S'-matrix elements. 

The structure of this paper is then as follows. In Sect. 2 we discuss the asymptotic 
fields and dynamics of QED in some detail. Then, in Sect. 3, we show how to correctly 
characterise charges and construct static charges. The solution for a charge in QED 
with an arbitrary, relativistic velocity is presented in Sect. 4. In Sect. 5 the field theo- 
retical implications of this approach are briefly reviewed and in Sect. 6 we present some 
conclusions and discuss the implications of this paper. 



2 Dynamics and Asymptotic Dynamics 

The starting point for our analysis of the physical, charged sector of QED is the familiar 
Lagrangian density 

£ = -\F„ U F» V + ^7 M (<^ - ieA^iP - mW> + \B 2 + Bd^ . (8) 
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Note that here we use Feynman gauge, B = —d^A^ 1 , for a discussion of other gauges, 
see, e.g., ||42||. There is a significant algebraic advantage to working in this gauge which 



leads to a much simpler account of our construction of charges. We insist, though, that 
our fundamental variables (which will be introduced in the next section) are themselves 
gauge invariant and hence insensitive to this choice of the Feynman gauge. In the 
companion paper II, we will indeed see the gauge invariance of our construction in 
many perturbative calculations. 

The extraction of the equal-time commutation relations from (§) is, at least, formally 
straightforward even in a fully interacting theory. What is much more difficult to do 
is to construct the general space-time commutators of these Heisenberg fields. For the 
B field, though, things are much better since in Feynman gauge it satisfies the free 
equations of motion (OB = 0). This fact, in conjunction with the known equal-time 



commutation relations, leads to the space-time commutators ||43|| : 

[B(x),B(y)} = 0, (9) 
[B(x),A li (y)]=idPD(x-y) (10) 

and 

[B(x)^(y)\ = eD(x-y)^(y). (11) 

where 

D(x -y) = -^-e(x° - y°)5((x - yf) (12) 

Z7T 

-S§?^ y)s ^ (x °- y ^' (13) 

and Uk = \k\. We stress that these commutators are valid for all times and, in particular, 
hold both at very early and at very late times. 

The significance of these commutation relations comes from the role of the B field 
in characterizing the physical states and observables of the theory. Given its trivial 
dynamics, we can expand this field in terms of its modes: 

B ^=l<^^M k) ^ +c[k) ^- (14) 

The condition on physical states is then that they are annihilated by the positive fre- 
quency partf] of B, i.e., 

c(A;)|phys) = 0, (15) 

for all k. Observables must preserve this condition, and hence must commute with c(k). 
From the space-time commutators ([TOD and ( [TT]) we find that 

[c(k),A,(y)] = -ik^ k -y (16) 

and 

[c(k)^(y)\ = -ie^(y). (17) 



5 The connection between this, by construction the positive frequency part of B, and the gauge and 
matter fields will be clarified in Eq. (37). 
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The first of these relations expresses the well known fact that not all components of the 
vector potential are physical: at the end of the day there should be just two physical 
photonic degrees of freedom per space-time point. In contrast, the second relation 
( |T7D is not usually interpreted as saying that the matter fields are unphysical. This 
is because this commutator contains an explicit dependence on the coupling constant. 
The widespread assumption then is that, in the remote past and future, the interaction 
vanishes and hence the asymptotic matter fields can be taken as physical. There is, 
though, no obvious reason why the large time limit of fll7D should vanish. Indeed, the 
long range nature of the electromagnetic interaction suggests that it is highly unlikely 
that a non-interacting regime is ever reached. In fact these arguments are usually 
expressed [(44J as an ad hoc assumption that the interaction is adiabatically switched off 
in the remote past and future. 

The appeal to such a 'switching off' mechanism for the interaction does not reflect 
the experimental set-up found in a typical scattering process. Rather, what happens 
is that the in-coming, or out-going, particles become widely separated and that it is 
this growing separation that effectively renders the coupling to vanish. In theories with 
no more than cubic interactions (such as spinorial QED) there is a simple argument 
first used (to the best of our knowledge) by Kulish and Faddeev to indicate that in 
theories with massless fields, this assumption is not true, i.e., the asymptotic dynamics 
in QED is not that of the free theory. 

The idea is very simple^: go into the interaction picture and take the large time limit 
of the interaction. If this limit is zero, then the asymptotic dynamics is that of a free 
theory. If the interaction does not vanish, then we simply cannot take the asymptotic 
dynamics to be free. 

In QED this works as follows. The interaction Hamiltonian in the interaction picture 
is given by 

H int (t) = -eJ d 3 x 4(t, x) J f "(t, x) , (18) 

where the free matter current is J M (t, x) = ^ (t, x)^ip(t, x) and we are using the su- 
perscript f to distinguish the free fields from the corresponding interacting or Heisenberg 
ones. To further fix our notation, we take as our free field expansions 

^ {x) = I ~m ^ s > s ^ e ' ip ' x + sy(p)j p - x ) , (19) 



where E p = y\p\ 2 + m 2 , a sum over s is understood and 

t rfih 1 

4(x) =jj^y 3 ^- k {^(k)e- lkx + a\{k)e^} . (20) 

We note here the basic space-time commutator for the free vector potential 
[Al(x),A j(y)} = -ig» v D{x - y) . (21) 

6 Although to make this argument more rigorous and applicable to a wider class of theories a bit 
more care is needed. In particular we stress that all such limits should be understood as weak limits 
between appropriate, normalisable states. The full construction with extensive applications will be 
reported elsewhere Q. 
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Inserting these free field expansions directly into (|18|) results in eight terms which 
can be grouped according to the positive and negative frequency components of the 
fields. Each of these pieces will have a time dependence of the form e 1 ^ 1 where ip 
involves sums and differences of energy terms. The argument used in H] (see also ITU 



and the discussion in supplement S4 of P6[ ) is that, as t tends to plus or minus infinity, 
only terms with if) tending to zero can survive and thus contribute to the asymptotic 
interaction. After performing the spatial integration, and using the resulting momentum 
delta function, one finds that in a theory describing massive charges only four terms of 
the form tp = ±(E p+ k — E p ± Uk) can survive in the large time limit. The requirement 
that E p+ k — E p ± Uk ~ can only be met in QED because the photon is massless, 
in which case it implies that Uk ~ 0, i.e., only the infra-red regime contributes to 



the asymptotic dynamics. From this observation it can be shown [p|,|45| that the full 

interacting Hamiltonian (|i~8l) does not vanish asymptotically!], but has in fact the same 
asymptotic limit as 

=~ef d 3 x A% x)J»{t, x) (22) 

with 

j^t, x ) = J0f w/ {p)63{x - tp/Ep) ■ (23) 

The operators in this current are only contained in the charge density 

p{p) = £(&t(p, s)b(p, s) - d\p, s)d(p, s)) (24) 



which implies that the asymptotic current satisfies the trivial space-time commutator 
relation 

mx),r as (y)] = 0. (25) 

As such, this asymptotic current can be interpreted as effectively the integral over all 
momenta of the current associated with a charged particle moving with velocity p/E p . 
This does not vanish as t — > oo. We thus see that the asymptotic dynamics of QED is 
not that of a free theory. 

It is important to note here that by 'asymptotic dynamics' we are referring to the 
dynamics in the neighbourhood of time-like and null boundaries of Minkowski space- 
time. At space-like infinity things are very different: the fields vanish and we impose 
the condition that the local gauge transformations become the identity there. This 
distinction is made clear in the Penrose diagram [[|7| of Figure |TJ 

What this argument has shown us is that, contrary to the statements found in the 
usual approaches to scattering in QED, the basic commutator ([TT|) really does imply 
that the matter field is not physical even in the asymptotic regime. The Lagrangian 
fermion in QED should not be associated with anything physical such as an electron! So 
what can we identify with physical matter? 

Given that we are dealing with a gauge theory with an unbroken symmetry, so that 
the reduction in the degrees of freedom must be in the gauge sector of the theory, it is 
clear E3I that we cannot now try to reduce the degrees of freedom in the matter sector. 



7 If the photon is given a small mass, as an infra-red regulator, then the asymptotic Hamiltonian 
does vanish. 
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Figure 1: Penrose diagram of Minkowski space-time: the dashed lines are the worldlines 
of massive particles while the short arrows depict incoming and outgoing radiation. We 
may only demand that fields vanish at space-like infinity, i Q . 

Rather, a procedure is needed for finding a gauge invariant combination of both matter 
and gauge fields — we call this general process dressing the matter fields. This will be 
the topic of the next section. Before moving on to that construction, though, it is useful 
to continue our previous analysis of the asymptotic form of the matter and gauge fields. 

It is important^ to realise that there are two clear steps in the extraction of the 
asymptotic form of the fields entering into our theory. 

• First we go from the Heisenberg fields to the interacting fields via the usual time 
ordered exponential of the interacting Hamiltonian: 

Oi(t,to) = Texp {-ij'dTH^r)^ K (t)f exp (ij'dr H int (r^j , (26) 

where T and T signify time and anti-time ordering respectively. This transforma- 
tion requires the introduction of an arbitrary time to at which the fields in the two 
pictures agree. This time is arbitrary at this stage because if we were to directly 
transform back to the Heisenberg fields all reference to to would vanish. At low- 
est order in the coupling we find that the relation between a Heisenberg operator 
On(t) and its interaction picture counterpart is 

Ox(t,to) = K (t)-i f dr[H int (r),0 H (t)}. (27) 

Jt 

8 In H and |t6j the form of the asymptotic fields were found by just translating from the interaction 
to the Heisenberg picture, this resulted in a rather ad hoc argument for why the contribution from the 
lower time to could be ignored in the Heisenberg fields. 
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The asymptotic dynamics for the Heisenberg fields is now found by using the 
asymptotic form of the interaction Hamiltonian to translate from the interaction 
to the Heisenberg picture. At lowest order in the coupling we get 

O^(t,t ) =Oi(t,to)+i rdr[iC(T),Oi(Mo)] 

Jto 

= H (t) +ifdr [H&(t) - H int (r), O u (t)} . (28) 

Jtn ' 



Note that there is a t dependence in this expression for the asymptotic fields at time t 
due to the potential contribution from the lower range of integration. This dependence 
is governed by the value of H^ t (t ) — H int (t ), so if we set £ = i 00 this will vanish by 
the construction of the asymptotic interaction. 

In practice, though, we do not want to start in the Heisenberg picture, rather we 
wish to go directly from the known fields in the interaction picture. What this argument 
then tells us is that we are free to ignore any explicit contributions from the overlap 
time to as long as this is taken to be large. 

The transformation from the interaction picture to the asymptotic Heisenberg picture 
is implemented by the usual time ordered product of the Hamiltonian (f2^). Using the 
commutator relations (|25|) and ([H|), this can be written as the product of two commuting 
terms: 

TexJifdrHZir)) = expU fdr H^(r)) expU fdr, [dr 2 [^(n) , ff£(r 2 )]) . (29) 

\J — oo / \J— oo / \ J —oo J —oo J 



From (|25| ) and fl2"ID, the double commutator [[H^ t (ri), H™ t (T 2 )], A 1 (x)] = 0, and hence 
the asymptotic vector potential is given by (t = x°) 

A-(aO = ex P (i£driy- t (r)) A^x) ex P (-i/W H&(rj) 

= 4(x)-e /* drd 3 yD(T-t,y-x)J^(r,y). (30) 

J —oo 

Note that it follows from this that 

= -eJ-(x) , (31) 

and hence we see that the asymptotic vector potential is made from a free part plus the 
field generated by the non-trivial asymptotic current, i.e., a response to the presence of 
any other charges. We also have from (|2T| ) and (p5|) that the asymptotic fields obey the 
same commutator as the free fields 

[A™{x),A™{y)] = -ig» v D{x-y). (32) 

The asymptotic form of the matter field follows in much the same way, but now both 
terms in (|29|) contribute. The resulting field is |§ 

^(x) = J ^D(p,t) {b(p,s)u s (p)e-^ + d\p,s)v s (p)e^} (33) 
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where the distortion operator D(p,t) factors into two commuting terms (reflecting the 
similar factorisation in Q25D): D(p,t) = D p haseip,t)D S oh{p,t)- The phase component 
-D P hase(p, t) does not depend on the gauge field (it is in fact gauge invariant) and only 
contributes to the phase of the operator (see the discussion in Supplement S4 of P5fl), 



as such we will neglect its contribution to the asymptotic field. We do note, though, 
that this phase depends on the presence of other charges in the system and hence enters 
into the distortion operator only at order e 2 and above in the coupling. 

The soft component of the distortion operator is given by 

^-^{^M s t^-'-t^)\ • (34) 

V soft / 

where we have made explicit that only the soft (infra-red) regime contributes to the 
asymptotic fields. Due to this distortion operator, the asymptotic field ip™ does not 
create or annihilate particle states. That is, if we were to define (as we would want to 
in the LSZ-reduction of the S-matrix) the operator 

b(q, s, t) := J d 3 x -J^'( q )^(x)ft** , (35) 



then rather than recovering the time independent, one particle annihilation operator, 
we get (neglecting the phase) 

b(q, s, t) = D soh (q, t)b(q, s) . (36) 

This operator does not have a single particle like interpretation, indeed it does not even 
act in a Fock space as it is a coherent state operator. The perturbative expression of 
this is the infra-red problem in QED: one cannot extract poles from the on-shell Green's 
functions of the theory and hence an S- matrix cannot be constructed ||TT|,|49|j. Such 
observations have led to the view among theorists || that no particle description is 
possible for the electron. Our point of view on this is that the matter field is simply not 
physical and we should not be shocked if its asymptotic limit is also unphysical: it is 
far too early to give up on the electron as a particle! 

In order to identify just what does describe the physics in the asymptotic regime 
we minimally need to find fields which commute with the B field. To this end we 
observe that at early and late times the B field can be written solely in terms of the 
free electro-magnetic potential: 

B( X ) -> = -0"4(s) . (37) 

This result follows from (|30|), and the conservation of the asymptotic current if 



hm f d 3 y ±-6(a - \y\)J?(a, y) = . (38) 

a 00 J \y\ 



To demonstrate this, we insert the explicit expression ( |23"| ) for the asymptotic current 
into (|38|) . This generates a product of delta functions whose arguments, it can be 
immediately seen, cannot both vanish for massive matter fields. 
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Hence the modes of the I?-field responsible for identifying the physical states (|15|) 
are asymptotically given by the modes of the free vector potential as 



c(k) = ik ■ a(k) . (39) 

Thus, the modes of the free fermionic field are physical in the asymptotic regime. 
Hence d^(q,s) |0) is indeed an asymptotic, charged one-particle state. But it is not 
the asymptotic limit of the Lagrangian matter field. Given that this charged state is 
created by a free creation operator, one might now ask if it indeed has an associated 
electromagnetic field. The important point to note is that the asymptotic potential (pOj) 
is not that of a free theory. In particular, using flT2|) and ([ft]), we see that 

dt(g, s)} = -± k / (g, 8 ) . (40) 

V (<Z ' x ) ~ Q x 



We recall from Chap. 14 of [50 that the classical electro- magnetic potential associated 



with a charge moving with momentum q is, at large distances, given by 

^class ( x ) = ^ . (41 ) 

" 47T J( q . x )2 _ q 2 x 2 

Thus the electromagnetic potential associated with the state cft(q,s)\Q) indeed corre- 
sponds, at large distances, to that of a particle moving with momentum q. Hence 
the asymptotically physical operator d>(q, s) creates both a charged particle state with 
momentum q and its associated electromagnetic potential. 

To summarise: we have seen that (f(q, s) is physical and creates the correct electro- 
magnetic field expected for a moving charge. It is, however, not directly related to the 
the asymptotic limit of the matter field ip. Is there then a (gauge invariant) operator 
which has cft(q, s) as its asymptotic limit? This is an important question since if there is 
such a field, then this would correspond asymptotically to an electron. The construction 
of such a field will be the subject of the next section. 



3 Static Charges 

We have seen that it is indeed possible to construct gauge invariant, charged particle 
states at asymptotic times in QED. In fact, these states are elements of the familiar 
Fock space associated with the modes of the free fermionic field. In this sense, QED is 
no different from any other theory with purely massive fields. Where QED differs from 
other (non-gauge) theories is in the fact that these fields are not the asymptotic limits of 
the matter fields from which the Lagrangian of the theory was constructed. The neglect 
of this simple fact generates the infra-red problem, prevents the construction of S'-matrix 
elements and leads to the abandonment of a particle description of charges. In order 
to motivate our general construction of the fields which do asymptotically correspond 
to the charged particle states, we will, in this section, reanalyse Dirac's construction of 
the static charged field and show in what sense it correctly describes a static charged 
particle. 
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Dirac's construction of the field ipn raises many fundamental questions that we need 
to address if we are to make progress in our aim to describe charged fields in both QED 
and, more generally, in the standard model. For example, it is far from clear how unique 
this construction is or how it should be extended to a non-abelian theory such as QCD 
where the chromo-electric field of a static charge is not known a priori. 

Before analysing the uniqueness of the construction, let us first of all see how a static 
interpretation of ipo emerges even when the fields are not infinitely massive. To this 
end we define the gauge invariant annihilation operator b(q, s, t, 0) by 

b(q, s, t, 0) := J d 3 x -^^{q) exp ( - *e^) (x)^(x)e^ . (42) 

The momentum q is taken to be on-shell in this definition. The additional label '0' in 
this operator reflects the conjectured static nature of Dirac's construction. 

At 0(e), the proposed annihilation operator ( [12]) is the sum of two terms, one from 
the matter as seen in Sect. 2 and one coming from the expansion of the dressing: 

b^(q,s,t,0) = Jd 3 x -l=u^{q)^{xy^-ie Jd 3 x -l= u ^( q )^( x )^( x ) e ^ , 

(43) 

where the superscript signifies that we only retain terms up to order e. Starting in the 
interaction picture, then transforming to the asymptotic Heisenberg picture will result in 
the first term in (f43l) becoming the 0(e) distorted annihilation operator fl36|) . Given that 
the second term in (|4"3|) is already at 0(e), its asymptotic limit in this approximation 
will simply be the large time limit of the operator expressed in terms of the free fields 
(0) and (f2T]). Using the identity 



rdH_ J_ fk_a^ _ ^m e ^\ f (44) 
J (2tt) 3 2u k V k 2 k 2 J 



and the Kulish-Faddeev argument discussed in Sect. 2, the second term in (53) can be 
shown to have the asymptotic form 

^ I # l t (^e-^ - ^>e-/*,) . (45) 

soft 

Combining the two contributions to b^(q, s, t, 0), we get, in the large time limit, the 
modified distorted annihilation operator 



&M(g,5,t,0)->&(g,s)jl+e J 



d 3 k 1 



soft 



(2vr) 3 2u k 



k 2 q-k ' ( 



k-a ] (k) _ q- aK k ) \ r itk- q /E a 
k 2 q-k 



(46) 



In general this distortion does not vanish and thus, even though we are dealing with a 
gauge invariant field, it does not allow for a charged particle interpretation. However, if 
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the 4- momentum q is at the static point in the mass-shell, i.e., when the four-momentum 
is mr]^, where r\ is the unit temporal vector, then the distortion terms are 

_ r _fk_ _L_ / fc-g(fc) itWk k-a^k) 
J (2tt) 3 2cu fc \ k 2 k 2 

soft 



(c(k)e- ituJk + c\k)e ituJk ) , (47) 



soft 



(2vr)3 2ul 



where we have used the asymptotic identification ( p9|) of the modes of the B field. As 
we have seen, these modes commute with the free annihilation operator b(q, s) and with 
themselves. Hence, using the definition of physical states, fll5"|) , we see that between 
physical states the distortion operator reduces to the identity operator at the static 
point in the mass shell. 

The above analysis shows that matter dressed a la Dirac, t/>£>, has a particle inter- 
pretation. This interpretation we stress only holds at the static point in the mass shell. 
Before discussing the exponentiation of this result, or indeed its extension to moving 
charges, we need to understand the uniqueness of the construction. To this end, we need 
to step back and ask where this solution came from, i.e., how should we characterise the 
dressing so that it describes a charged particle with a well defined velocity? 

Let us initially look at the simpler situation where we have a heavy matter field, 
ip(x). This field can be thought of as the infinite mass limit of either a scalar or Dirac 
field which is, for the moment, not coupled to a gauge theory. In this limit velocity is 
superselected J51] and the equation of motion for the field has the universal form 



u ■ d(p(x) = , (48) 

where u is the associated 4-velocity of the heavy particle. This equation is simply the 
statement that the field is constant along the world line of a particle moving with 4- 
velocity u. Indeed, if we parameterise the world line of a particle moving with 4-velocity 
u ix = ry(jj _|_ (-where r\ is the unit time-like vector introduced above, v = (0, v) is the 
velocity and 7 = (1 — |t>| 2 ) -1 / 2 ) as 



x 



(s)=x* + (s-x )(7i + v) ,i , (49) 



then (13) implies that for arbitrary s, 

<p(x(s)) = <p(x(0)) . (50) 



If the heavy matter field is now minimally coupled then the equation of motion ([48|) 
becomes 

u-D<p(x) = 0, (51) 

where = — ieA^. There is now no heavy charged particle interpretation to this 
equation. However the field ip should not, and indeed cannot, be identified as a physical 
field since it is not gauge invariant. The lack of a simple particle interpretation to the 
field ip is thus not a serious problem since what we need to do first is to dress this field. 
We will demand that the dressing is such that, in addition to restoring gauge invariance, 
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it should also ensure that a (heavy) particle interpretation for the dressed matter field 
holds. 



We define the gauge invariant charged matter field by 

$(x) = hT 1 (x)<p(x), (52) 

where the dressing h~ 1 (x) is such that under a gauge transformation we have 

h-\x) -> h-\x)e- ieeix) . (53) 

That is, a physical particle corresponds to the heavy matter field dressed with some 
'brown muck', whose exact form we will now clarify. For this field to correspond to a 
heavy charged particle we must have u ■ <9$ = 0. This follows from (|5lD if the dressing 
satisfies 

u-dh- 1 {x) = -ieh- 1 {x)u- A{x). (54) 



We call this equation the dressing equation. Equations (|53| ) and ( p4j ) are the fundamental 
requirements on any description of charged particles and are central to what follows. 

Although our motivation for this equation emerged from an analysis of the heavy 
matter sector, we demand that, more generally, it holds for the dressing of any field that 
is asymptotically corresponding to a charged particle moving with velocity u. There 
are two arguments for this: it is well known that the asymptotic dynamics of QED is 
governed by soft photons for whom any electron is heavy; secondly, it can be shown |10| 
that the asymptotic interaction Hamiltonian vanishes for the propagator of the dressed 
fields which satisfy ([54]) if one is at the right point in the mass-shell. We now specialise 
to the static situation and leave the general case to the next section. 



In the static situation the dressing equation fl54j) becomes 

d Q hr x {x) = -ieh' 1 (x)A (x) . (55) 

The first thing to note about this is that Dirac's proposal (||) does not satisfy this 
equation! Indeed, using the identity 

d„e° = e°(d,0 + \[d,0 ) 0}) (56) 

where O is an arbitrary operator whose commutator [d^O, O] is a c- number, we find 
that 



. diAi. \ ( . diAi. \ dodjAj e r d A k 1 \ 

- - vr(s) J = — ex P { ~ "vrW) (^-(*) + 2 J (27p ul) ' 

and even if we ignore the field independent term in the last factor, we do not solve (|5 

To understand the relationship between Dirac's dressing and the static version of 
the dressing equation we need to solve equation fl55|). As is well known, the solution to 
equations like (|55"D have the form of an anti-time ordered exponential: 

/i _1 (x, a) w Texp ^ — ie J A (s,x)ds S j, (58) 
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where a is, as yet, undetermined. Although this clearly solves (p5|), it does not satisfy 
the gauge transformation property Q53D which is essential for a dressing. Indeed, under 
a gauge transformation we have 

h-\x, a) -> e iee ^h~\x, a)e" iee(x) . (59) 

There is no choice for the parameter a for which 9(a, x) vanishes (the choice a = ±00 
evaluates 9 at time-like infinity where there are no restrictions on the gauge transfor- 
mations, see Fig. [I]). 

We can compensate for this bad behaviour under gauge transformations, while still 
solving the dressing equation fl5"5]), by taking instead 

h~ 1 (x, a) = exp ^ — ie-^ L (a, x)*jT exp ( — ie J A (s,x)ds\. (60) 
This we can write as 

h~ 1 (x, a) — Texp ( — ie J {A (s, x) - (s, a;)| ds^j exp ( - ie^^-(x)^ , (61) 



where we have combined the two exponentials of (60) under one exponential, then 
written one part as a total derivative. We have neglected, for the moment, possible 
commutator terms which will be discussed in Sect. 4. The first term in this expression 
is now gauge invariant, a fact which becomes manifest when we write this solution as 

h^(x, a) = Texp (^ie J - > ^ t ° ( s j x ) ds^j exp ( — ie^^(x)J . (62) 

Thus we see a factorisation of the static dressing into a minimal part which is essential 
for gauge invariance (and is just Dirac's original proposal for the static dressing), plus 
an additional part which is separately gauge invariant. 

This general solution to the static dressing equation takes on a much simpler form 
in the asymptotic regime where a charged particle picture should emerge. Using the 
asymptotic space-time commutators fl32[) we see that 



= 0- (63) 

This commutator implies two important results. Firstly, the additional part of the 
dressing will commute with the electric and magnetic field operators, and thus the 
electromagnetic field associated with the static dressing is the same as that produced 
by the minimal, Dirac component. Hence we see that Dirac's argument for the form of 
the dressing was not sensitive enough to detect the additional term. Secondly, we can 
use this commutator to dispense with the anti-time ordering in the dressing (|62"D so that 
in the asymptotic regime 

h^(x, a) — ► exp [ie J (s, x) dsj exp ( — ie 1 * (x)j . (64) 

Using Gauss' law (@) in the asymptotic domain, and the identification (|30|), we see that 
the static dressing factorises in the asymptotic regime into the product of two terms. 
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One can be made out of the various matter contributions and is gauge invariant, it is 
analogous to the phase term in the distortion operator for the matter field. The other 
part of the dressing is constructed out of the free vector potential and, for reasons that 
will become immediately apparent, we call it the soft part of the dressing: 

Kk( x ) = ex P ( ~ ie ^r( x )) ■ ( 65 ) 



This is the most important component in QED and, as was shown in J4C§, its extension 
to QCD is the dominant glue around quarks at short distances. 

Armed with this result, we can now study the exponentiation of the soft contributions 
to the asymptotic annihilation operator Repeating the Kulish-Faddeev argument 
for the soft components of both the matter and dressing we find that 

b(q, s, t, 0) -> /i~ft(g, t)D soh (q, t)b(q, s) , (66) 

where 




soft 



(2tt) 3 2iu k 



which is the exponentiation of our previous result fl45|) . The combined effect of these 
distortion factors can easily be evaluated using the canonical commutation relations for 
the modes a^(fc) and the Baker- Campbell-Hausdorff formula. One finds 



^soft^^Arftte,^ = exp 



d 3 k 



(2tt) 3 2oj, 



k-a{k) _ g-a(k) \ itk . q/Eq 
k 2 q-k 1 



' soft 

f k-a\k) _ q-a ] (k) \ ltk . q/Ea 
k 2 q-k U 



(68) 



Note that the commutator in the BCH- formula vanishes here. Just as we saw in (47Q, 
this becomes the identity operator between physical states at the correct point in the 
mass-shell where the momentum q^ is at the static point in the mass shell: q^ = mrf '. 

To summarise: we have seen in this section how to characterise the dressing appro- 
priate to a charged particle moving with a given velocity. In the static case we solved 
the dressing equation and saw how Dirac's proposal for a static charge emerged as the 
minimal part of the static dressing. Using this construction we showed how the static 
annihilation operator emerges as the asymptotic limit of the statically dressed matter 
field. Having seen how to describe static charges, we now proceed to an analysis of 
moving charges. 



4 Moving Charges 

In this section we will show that the dressing needed to describe a charged particle 
moving with four-velocity is 

h-\x) =e- ieK{x) e- iex{x) , (69) 
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where the minimal part of the dressing is given by 

X(x) = §^ > (70) 

with Q 11 = (r) + vY{rj — v) ■ d — and the additional (gauge invariant) part of the 
dressing is 

K(x) = [ X \v + vr^f(x(s)) ds . (71) 
J±oo y ■ a 

In this last expression the integral is along the world-line of a massive particle with 
four- velocity parameterised as in fl49|) . Before deriving this result from the dressing 
equation ([54]), we will first verify that this dressing generates the correct electric and 
magnetic fields and demonstrate that its modes have the asymptotic limit appropriate 
to a charged particle moving with velocity v. 



In the asymptotic regime we can again argue from (^3|) that the additional component 
of the dressing does not affect the electromagnetic configuration. This allows us to define 
the minimally dressed charged field ip v by 

if> v ( x ) = e-^VW • (72) 

By construction, this is gauge invariant. It is also straightforward to see that this field 
is physical in the context of the -B-field formalism discussed in Sect. 2. In particular, 
from ([Uj) and (|TTD, it is easy to see that 

[c(k),ip v (x)) = 0. (73) 

In the static limit this dressed field reduces to Dirac's expression (||), i.e., ipo = ipD- 

Repeating the argument given by Dirac leading up to the expression (|6|) for the 
electric field of the static charge, the electromagnetic fields generated by ( |7^ ) are given 
by the equal-time commutators — ie[Ei(x), x(y)] et and —ie[B i (x),x(y)] et - In order to 
evaluate these we first need to make clear how 1/Q • d is defined in (^) . 

We take 

<^/(a0 := / d Z z G(x - z)f(z) , (74) 

where 

g-dG{x-z) =5 3 (x-z) (75) 

and we note from its definition that Q ■ d = V 2 — (v ■ d) 2 . Taking the Fourier transform 
of (|75|) we see that 

G <*>=-/SsMbr**- (76) 

This integral is computed in the standard fashion by first diagonalising the quadratic 
form in the denominator. The final expression for G(x) is then 

G{x) = ~ . 7 (77) 



Note that in the static limit this reproduces 
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The electric and magnetic fields generated by this dressing are now straightforward 
to calculate. For example, applying the operator ( |T2"D to a state |0) with no electric 
field, we find 

E^MvM = -~ ^ {x ~ y ) t^UvM • (T8) 

4?r - ?/| 2 + 7 2 (-u • (x-y)) 2 Y 

This is precisely the electric field associated with a charged particle moving with veloc- 
ity v. 

As we saw with the static charge, the interpretation of the operator ip v (x) as that 
creating a moving charge at the point x only holds in the infinite mass limit. More gen- 
erally then, we extend the definition of the static, gauge invariant annihilation operator 
(fUD to the moving case by defining the mode 

b(q, s, t,v) = J d 3 x -^u^iq) exp ( - ze|-|) {x)^{x)^ x . (79) 

We will now demonstrate that this is a particle annihilation operator at the right point 
in the mass-shell. 

At 0(e), b^(q, s,t,v) is the sum of two terms in which the free field expansions for 
the potential in the dressing can be used. Using the expression 

« ' * , x) = i (** JL f ^e"- - , (80) 



G ■ d v ' J (2tt) 3 2oo k V V ■ k V ■ k 

where we have introduced the notation that 

V = (77 + vY{r) - v) ■ k - jfc" , (81) 
we see that in the large time limit 



'V-a(k) q-a(k)\ itk 



q ■ k 



c 



g/E. 



V-a){k) q-a\k) 



e itk-q/E q 



V ■ k q ■ k 

In the distortion term we note that we can write (recall that k is on-shell) 

(f { J] + v Y(- n -v)-k-k il <f 

V ■ k q ■ k (77 + v) • k(r) — v) - k q ■ k 

(77 + vY g M k^ 
(j] + v) ■ k q ■ k V ■ k 



(82) 



(83) 



Hence, at the point in the mass-shell where q^ = m / ~f(r] + v ) M , the distortion operator 
becomes the trivial operator 



f d " k 1 ( k-a(k) itu)h k ■ gi(k) itU) \ 
J (2tt) 3 2u k \ V-k V ■ k ) 

d 3 k i 



(2tt) 3 2u k V ■ k 
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(c(k)e- ituJk +c j (k)e ituJk ) , (84) 



which vanishes between physical states. 



The exponentiation of this result for the soft contributions to both the distortion 
operator and the dressing, follows in much the same way as it did for the static case. 
Indeed, at large times we get 

b(q, s, t, v) -> h soft (q, t, v)D soit (q, t)b(q, s) , (85) 

where D so{t (q,t) is given in (0) and 

^ Av) .^U^(^^^-y^^)\ . ,8a, 

vsoft ) 

The combined effect of these distortions yields the overall soft distortion of the particle 
mode as 



v _ . * / f d 3 k 1 \{V-a(k) q-a(k)\ 

h so{t (q,t,v)D so{t (q,t) = exp^eJ^—^-y^--^) 

soft 



'V-a(k) q-a(k)\ _ itk 



e -*' E * 



V -flt(fc) _ q- o t (fc) 
V ■ k q ■ k 



itk-q/E, 



(87) 



Just as we saw above, this becomes the identity operator between physical states at that 
point in the mass-shell where g M = 7717(77 + v)^. 

This analysis has shown that the dressed matter ip v can be interpreted as a charged 
particle moving with velocity v at the appropriate point in the mass shell. In this sense, 
we see that a relativistic description of charged particle is indeed possible. This dressing 
fl69D for a moving charge can, of course, be derived from the dressing equation (0). Full 
details of this are given in the appendix. 



5 From Dressed Matter to the /S-Matrix 

In this section we want to show how the S'-matrix elements of QED can be extracted 
in terms of Green's functions of dressed matter fields. We will see that the usual per- 
turbative techniques can be applied. It should become clear that a correct application 
of the usual LSZ formalism (for a discussion with fermionic fields see, e.g., Chap. 5 of 
Ref. |5"2fl ) requires the use of dressed matter. 

Recall that in the usual LSZ route to the S'-matrix it is assumed that there is 
an overlap at large times between the interacting fields and the free fields. This is 
schematically illustrated in Fig. |2|. Since we have seen that this overlap does not take 
place between the usual matter fields in the QED Lagrangian, it is clear that we should 
expect problems in implementing the LSZ procedure and, indeed, one finds that the 
usual Green's functions do not have a good pole structure. The responses to this range 
from using coherent states and abandoning the pole structure |k|) to giving up the idea 
of S'-matrix elements and only considering inclusive cross-sections including all possible 
final state photons. We now want to show how the use of dressed fields in QED lets us 
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Figure 2: Sketch of a scattering process where the coupling is assumed to vanish for 
large times. 

follow the LSZ path to the S-matrix. The companion paper, II, will show in perturbative 
studies that good pole structures are obtained with the use of these fields. 

Consider the simplest case: one particle goes to one particle. The equivalent S- 
matrix element may be expressed as 

(out | in) = (out \tf(q)\ 0) , (88) 

where |0) is the vacuum state and (b) are just the traditional one-particle Fock creation 
(annihilation) operators. In the free theory we may express the creation operator by 

&t(g) = Jd 3 x 7^r^(x)u(q)e- tq - x . (89) 

Now it is normally assumed that for large times the usual Heisenberg fields approach 
the free fields. However, we have seen in Sect. 2 that this is not the case in unbroken 
gauge theories like QED. Rather, as was demonstrated in Sect.'s 3 and 4, we have 
at each point, q, on the mass shell an appropriately dressed field, ip v , such that in 
( |89| ) we may, if we neglect the unobservable phase effects, (weakly) replace ip v — ► ipf for 
t — > ±oo. (Note that for simplicity we do not explicitly insert renormalisation constants. 
Renormalisation will be studied in the explicit calculations of II.) 

We thus see that in QED we obtain 

<out|6 t (g)|0> = lim fd 3 x -^-(out\^l(x)\0)u(q)e- iq - x . (90) 
t-*-oo J 2E q 

In the standard fashion we may introduce an integration over time to obtain 

(out|fet( g )|0) = -Jd 4 Xj t ^(out\^ v (x)\0Mq)e-^ 

+ lim / d 3 x ^(out\i)l(x)\0)u(q)e- tq - x . (91) 
t^oo J 2E q 

As is usual this second, reduced term corresponds to disconnected graphs and may be 
dropped. 

The time derivative acts upon both the matter field and the exponential. Writing 
ipl = tpv'jo and using that (^ — m)u(q) = 0, we find 

(outl&tfaOlO) = i I d 4 a;(out|^(a;)|O)(0 + m)^—u(q)e- iq - x . (92) 

J %Eq 
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We note that this is the standard result, but the matter field is replaced by the appropriate 
dressed matter. 

We now want to continue by using an annihilation operator to replace the out state 
by the vacuum, using the counterpart of (jS9f), and taking time-ordering into account, 
we rapidly find 

(out|in) = (0\b(p)tf(q)\0) 

= Jd A x Jd 4 y^ y - m)(0\T (Mv)Mx)) \0)(i^ x +m) 

x^r«(g)e-^^r«(p)e^, (93) 

where p (and v') correspond to the momentum (velocity) of the outgoing charge. Clearly 
this can be generalised to more interesting S'-matrix elements. 

We now want to reexpress ([53) in terms of a Gell-Mann-Low formula. The usual 
method relies upon the existence of a unitary operation linking the Heisenberg and the 
free fields. We have already seen in this paper that the free matter fields are not unitarily 
related to the equivalent Heisenberg fields, but for the dressed fields we indeed have 

Mt,x)=U- 1 (t)tf,{t,x)U(t), (94) 

where U is the usual operator, made up from exponentiating the interaction Hamiltonian 
expressed in terms of free fields. We stress that this relation only holds at the right point 
on the mass shell appropriate to the dressing. 

Once we have this relation it is straightforward to repeat the usual procedure to 
obtain the standard Gell-Mann-Low formula but with the charged matter fields replaced 
everywhere by appropriately dressed matter. With this, we are in a position to perform 
perturbative tests. These are presented in the companion paper, II. 



6 Conclusions 

In this paper we have shown how to describe charged particles in relativistic QED 
through a process of dressing the matter with the appropriate electromagnetic cloud. 
The essential ingredients which we used to obtain the precise structure of the dressing 
were gauge invariance (|53|) and the kinematical requirement which we call the dressing 
equation ([54]). 

We solved these two conditions and obtained a structured dressing composed of two 
factors: a minimal component which had the correct gauge transformation properties 
and an additional, gauge invariant part which was necessary to fulfill the dressing equa- 
tion. In this paper we have focused primarily on the minimal part of the dressing and, in 
particular, we have shown that it removes the soft distortion factor which has been pre- 
viously claimed to prevent a particle interpretation of QED. This has been demonstrated 
for a charged particle moving with an arbitrary velocity. 

The vanishing of the distortion factor shows that the long range interactions of 
QED, due to the masslessness of the photon, are encoded in the dressing rather than 
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any residual interaction term in the Hamiltonian. An important consequence of this is 
that the asymptotic Hamiltonian becomes the free Hamiltonian for our dressed fields. 
This means that the LSZ-formalism can be directly carried out without any use of a 
dubious adiabatic 'switching off' of the QED coupling. 

The process of dressing a gauge dependent variable so as to make it gauge invariant is 
a very general procedure. How one would want to dress the fields that form a bound state 
might be very different from how individual charges should be dressed. In this work we 
are only interested in dressing a charged field so that it corresponds to a charged particle 
moving with the appropriate velocity at either early or late times. The interesting topic 
of how to describe, say, a bound state such as positronium, that is either entering or 
leaving a scattering process, we leave to future work. 

Turning to QCD, where we want to construct colour charged quarks and gluons, 
we note first that our two fundamental requirements on the dressing have immediate 
non-abelian extensions. A simple algorithm for constructing the minimal, soft gluonic 
dressing around a quark can be found in pi|| . The potential between two so-dressed 
quarks was studied in [[5J. There it was shown that this dressing generates the anti- 
screening interaction of QCD which is responsible for asymptotic freedom. We have thus 
identified the non-abelian extension of the soft dressing as the dominant glue around 
quarks. A corollary of this result is that we see that the anti-screening interaction takes 
place between two separately gauge invariant constituent quarks and there is as yet no 
sign of any flux tube formation. 

In the companion paper II we will take the results of this paper and submit them to 
a battery of tests. There we will see that the particle of our dressed fields also emerges 
from perturbative calculations. The Green's functions of these fields will be shown to be 
free of on-shell infra-red divergences and to have the pole structure expected of physical 
particles. Detailed studies of the ultra-violet behaviour of these fields and their Green's 
functions will be presented and it will be shown that the renormalisation programme 
can be carried out without difficulties. 
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A Appendix 

In this appendix we show how the dressing for a moving charge (|69|) follows directly 
from the gauge transformation properties (|3|) of the dressing and the dressing equation 

©■ 

Before presenting the solution, we need to briefly reexamine Dirac's general for- 
mula (§) for a gauge invariant field. In line with @, we take as an ansatz for the 
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electro-magnetic cloud an exponential of the form exp(— iex{x)) where, under a gauge 
transformation 

X (x) -> X {x) + 9{x) . (95) 

This is precisely the type of transformation that Dirac investigated in (|3|) and we would 
be tempted to accept as the general solution 

X (x)= (d A zr{x-z)A,{z) 1 (96) 



where f^{x — z) satisfies <9 M / M (x — z) = 5 A (x — z). But this only implies ( [9"5] ) if no 
surface terms arise when we integrate by parts. The restriction on the local gauge 
transformations to those that vanish at spatial infinity is quite natural as finite energy 
restrictions impose a 1/r fall-off on the potential. However, as discussed around Fig. |], 
no such restriction applies to the fields at temporal infinity. Thus, to maintain gauge 
invariance, the form of f^{x — z) must be such that no surface terms arise at large times. 

As it stands, we can only infer from this that f°(x — z) should be zero outside of some 
bounded region in the z°-direction. To get more from this, we follow Dirac's proposal 
for the static charge and take 

X (x) = J d A z G(x - z)g"Ap{z) , (97) 

where is a first order, differential operator and Q ■ dG{x — z) = 5 4 (x — z). In order 
to avoid the surface terms that would obstruct the gauge transformation properties of 
the dressing, we must have that the operator Q ■ d cannot involve any time derivatives. 
Given this restriction, we see that G(x — z) = 5(x° — z°)G(x — z) where 

g-dG(x-z) = 5 3 (x-z). (98) 

We shall, for convenience, write x as 

X (x) = Jd 3 zG(x-z)g»A,(x ,z):=^(x) 1 (99) 

which mirrors the notation adopted in the body of this paper (cf. Sect. 4). So far, of 
course, Q remains undefined here. 

We now make a more general ansatz, based upon our experience in the static case, 
i.e., that the dressing factorises into the product of two terms: 

h~\x) =e- ieK(x) e- iex{x) , (100) 

where K is gauge invariant and x ~~ * X + @ under the gauge transformation (|J). The 
second term is thus a minimal part of the dressing, essential for the gauge transforma- 



tion properties of the dressing. The first term in ( 100 ) is then an additional part of 
the dressing. Here we write this as a simple exponential, rather than a path-ordered 
exponential as we did for the static example in (|58f) , since we are taking our fields to be 
the asymptotic fields discussed in Sect. 2, and the commutators of these fields will be 
seen to allow this ansatz. We note that all of the fields in the rest of this appendix are 
asymptotic fields and we omit the 'as' superscript in what follows. 
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For the additional part of the dressing it is enough in this abelian theory to also take 
K to be linear in the fields and we find that it is sufficient to make the expansion 



eK{x) = eK^x) + e 2 K 2 {x) . (101) 

The dressing equation ( 54]) can then be expanded in powers of the coupling to give the 
two equations: 

(77 + vYd^K, + X ) = (V + v Y A n , (102) 

and 

i(rt + vYd,K 2 = \{ v + vY(d,{ X , K x ] - [Ap, K x + x \) • (103) 
The solution to ( |102| ) is 

K x {x) + X (x) = [ X (V + vYAp(x(s)) ds + X (x(a)) , (104) 

J a 

where the integral is along the world line parameterised by ( [49]) with starting point a. 
The constant term %(x(a)) is needed to ensure that K is gauge invariant. This means 
that we can write K x in the manifestly invariant form: 

K x {x) = £ (r l + vr(Ap(x(s))-d,x(x(s)))ds 

\v + vY^f(x(s))ds. (105) 

The lower range of integration found in K\ should, in this asymptotic regime, have 
no physical significance. This will be ensured if the derivative of the additional part of 
the dressing with respect to a is zero between physical states. This would be the case 
if this derivative either vanishes or is constructed from the I?-field. From the form of 
(|105|) we see that it will not vanish. Therefore we must show that it is made from the 



.B-field. We thus postulate that it can be written as 

^K(x) = (v + vr(^(x(a)) + e-^(x(a))), 

= (r) + vr^(x(a)). (106) 

We now investigate what this implies for K\ and K 2 . 
From equations ( |105| ) and (|106|) we see that 

= (v + vr^f(x(a)). (107) 

Hence, from the anti-symmetry of F Ufl and the fact that Q u is a first order differential 
operator constructed out of d v , rf and v u , we must have that 

Q v = -d u + {r] + v) u {ar]-d + (3vd). (108) 
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The constants a and (3 can be fixed by the requirement that Q ■ d contains no time 
derivatives. It is easy to see that 

g ■ d = (a - 1)(t7 • df + {a + /3) (77 • d)(v ■ d) + (3(v ■ df + V 2 . (109) 

Hence we must have that a = 1 and (3 = — 1. From this we see that 

g v = (rj + v y(T)-v) -d-d u . (110) 

Given this result we obtain 

Ki(x) = £\ V + vr^f(x(s))ds. (Ill) 
Using the commutators ( |63D it follows that 

[K 1 (x),jj±(x)] = Q (112) 



and hence the a-independence of the construction will follow from ( |106|) if 

^) = (^^Wa))- (H3) 

This equation and (|103|) determine the form of K2. 

Given that the (asymptotic) current, J, is gauge invariant, we can immediately solve 
(|TT3D to get 

/a J 
(rj + v) v — v (x(s)) ds + a-independent terms. (H4) 
00 y • o 

Putting this expression into ( TO we see that the a-dependence cancels on the left-hand 
side. This means it cannot survive the commutators in ( |103| ). This cancellation is far 
from obvious. To show that it does happen we note that 

d 3 k 1 
(2tt) 3 uj 2 -{k-v) 2 

<! " k 1 lfc ^°- a )cosK(x°-a)) (115) 



( V + V y[A IM (x),K 1 (x)+x(x)]=ir 



(2vr) 3 ul-(k-vf 

_fk z^v pikMx o. a) gin ( ^ fc(x o_ a)) 



(2tt) 3 u k {ul-{k-v) 



and 



[xW, k M] = - [is f*!L _ R L_ e .-(^-»» ( cos Ml »_ s)) 

+ sin (io k (x°-s))){llQ) 



Hence 



(77 + V ydMx),K x (x)\ = d \x{x), K x {x)\ (117) 

— o ^(* -a) cos (^( x o_ a )) (118) 



(2tt) 3 uj 2 k -{k- v) 2 

-ilwT* i^'il ^ e^°^sin(^(x°-a)) 
(2tt) 3 u) k {u){-{k ■ v) 2 ) 
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which implies that that a-dependence has indeed been killed off. Thus from ( |103| ) we 
see that 

r d 3 k 1 

,(„ + „r Wx) = -In • a") 

Combining this equation with ( |113| ) we can solve for K 2 to get 



k 2 {x) = -§7-^ • x/ 7^ 2 fu — 9 + 1 (77 + ^r^o^)) rfs • ( 12 °) 



This together with (|111|) , ( |99"|) and (|1 10|) is the general solution to the dressing equation 
and the requirement of gauge invariance. Since, by construction, this is independent of 
a, we can set a = —00 and hence obtain the result promised at the start of Sect. 4. 
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